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ABSTRACT 

Based  upon  the  ray  theory/  we  develop  a  systematic  method  to  obtain  an 
equation  of  K-dV  type  with  variable  coefficients  for  the  evolution  of  water 
waves  in  a  channel  of  nonuniform  cross  section.  Examples  for  channels  with  a 
nonuniform  rectangular  and  triangular  cross  section  are  given.  The  fission  of 
solitons  in  a  triangular  channel  with  a  shoal  is  studied  by  the  inverse 
scattering  method  and  also  numerically.  A  general  Green* s  law  for  the  decay 
of  wave  amplitude  in  a  channel  with  arbitrary  cross  section  is  derived. 
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SIGNIFICANCE  AND  EXPLANATION 


In  this  paper  we  propose  a  new  equation  for  the  study  of  nonlinear  water 
waves  in  a  channel  of  variable  cross  section.  We  also  use  both  analytical  and 
numerical  methods  to  investigate  the  problem  of  a  solitary  wave  climbing  up  a 
shoal  in  the  triangular  channel.  A  solitary  wave  is  a  wave  with  a  single  hump 
moving  with  a  constant  velocity.  It  is  shown  that  under  certain  conditions 
the  solitary  wave  may  split  into  several  similar  waves.  Furthermore,  a 
general  law  to  describe  the  attenuation  of  a  wave  in  a  channel  is  also 
derived. 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MRC,  and  not  with  the  authors  of  this  report. 
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WATER  WAVES  IN  A  CHANNEL  OF  VARIABLE  CROSS  SECTION 

*  t 

Xi-Chang  Zhong  and  M.  C.  Shen 

1.  Introduction. 

The  K-dV  equation  originally  derived  for  water  waves  over  a  uniform 
bottom  has  been  studied  extensively  since  the  inverse  scattering  method  was 
discovered  by  Gardner,  Greene,  Kruskal  and  Miura  (1974)  for  the  solution  of 
the  equation.  Recently  there  has  been  growing  interest  in  K-dV  equations  with 
variable  coefficients,  which  appear  in  the  study  of  water  waves  over  a  bottom 
of  nonuniform  depth.  The  K-dV  equation  for  a  two-dimensional  variable  bottom 
was  derived  by  Kakutani  (1971)  and  Johnson  (1973),  and  that  for  a  rectangular 
channel  with  variable  width  and  depth,  by  Shuto  (1974).  A  review  of  the 
recent  developments  regarding  K-dV  equations  with  variable  coefficients  may  be 
found  in  an  article  by  Miles  (1980).  So  far  all  the  work  on  K-dV  equations 
for  water  waves  over  a  variable  bottom,  one  way  or  another,  is  related  to  the 
equations  for  a  rectangular  channel,  and  the  method  of  derivation  cannot  be 
extended  to  channels  with  arbitrary  nonuniform  cross  section.  On  the  other 
hand,  the  equations  derived  do  not  depend  upon  time  explicitly,  and  conse¬ 
quently  an  initial-value  problem  cannot  be  posed.  Nevertheless,  some  inter¬ 
esting  phenomena  have  been  observed  from  a  study  of  these  equations  concerning 
the  effect  of  a  perturbation  on  a  solitary  wave  propagating  in  a  channel  of 
variable  depth.  A  single  soliton  may  split  into  a  finite  number  of  solitons 
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if  the  depth  of  the  channel  decreases  from  one  constant  depth  to  another 
(Tappert  and  Zabusky,  1971;  Johnson ,  1973).  Furthermore ,  the  generation  of  a 
shelf  is  observed  behind  a  solitary  wave  as  it  propagates  in  a  rectangular 


channel  with  a  slowly  varying  bottom  as  a  perturbation  about  a  constant 
depth  (Ko  and  Kuell,  1978;  Miles,  1979;  Knickkerbocker  and  Newell,  1980). 

The  purpose  of  this  paper  is  to  propose  a  general  K-dV  equation,  which 
applies  to  any  channel  with  variable  cross  section.  The  method  of  derivation 
used  here  is  essentially  a  specialization  of  the  procedure  developed  by  Shen 
and  Keller  (1973).  Their  work  was  motivated  by  the  results  due  to  Choquet- 
Bruhat  (1969)  on  ray  method  for  nonlinear  partial  differential  equations  and 
by  Keller's  linear  ray  expansion  (1958)  for  water  waves  over  a  variable 
bottom.  In  terms  of  the  time  variable  along  a  ray,  the  K-dV  equation  is  now 
time- dependent  and  can  be  used  to  study  initial  value  problems.  Furthermore, 
there  is  great  flexibility  at  our  disposal  to  choose  the  appropriate 
independent  variables  in  the  K-dV  equation,  and  a  general  Green's  law  for  the 
decay  of  wave  amplitude  in  a  channel  of  variable  cross  section  is  also 
derived.  Based  upon  the  K-dV  equation,  we  shall  use  inverse  scattering  method 
to  study  the  fission  of  solitons  as  a  solitary  wave  climbs  up  a  shoal  in  a 
triangular  channel.  The  results  are  then  confirmed  by  numerical  methods. 

In  Section  2  we  formulate  the  problem  and  derive  the  general  K-dV  equa¬ 
tion  by  means  of  the  ray  method.  In  Section  3  we  specialize  our  result  to  the 
equations  for  rectangular  and  triangular  channels.  The  analytical  and  numer¬ 
ical  results  for  the  fission  of  solitons  in  a  triangular  channel  are  given  in 

l  Section  4,  where  the  change  of  the  channel  width  or  depth  is  considered.  In 

I 

I  Section  5,  a  discussion  of  the  results  is  given  and  the  general  Green's  law  is 

derived. 
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2.  Derivation  of  the  K-dV  equation. 

We  consider  the  motion  of  an  inviscid,  incompressible  fluid  of  constant 

*  *  *  * 

density  under  gravity  in  a  channel  with  a  boundary  defined  by  h  (x  ,y  ,z  )  = 


0,  where  z  is  positive  upward  (Figure  1),  The  governing  equations  are 

*  *  * 


u  * 

+  v#+w#  =  °  , 

(D 

X 

y  z 

* 

*  *  *  *  *  * 

* 

u  * 

+  uu#  +  vu*  +  wu* 

=  -p  */p  • 

(2) 
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X  y  y 
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* 

*  *  *  *  *  * 
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* 

w  * 

+  uw#  +  vw#+ww^ 

=  -p  */P-g, 
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x  y  z 
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subject  to  the  boundary  conditions 

#  *  * 

*  *  *  -\ 

*  *  +  u  n  # 

+vn#-w  =o  j 

*  *  *  *  * 

(5) 

t  X 

y  > 

at  z  =  n  (x  ,y  ,t  ), 

p  -  o  J 

(6) 

*  * 

*  *  *  * 

* 

u  h 

#  +  vh#+wh#  =  0 

at  h  =0. 

(7) 

x  y  z 

★  *  *  it 

Here  (u  ,v  ,w  )  is  the  velocity,  t  is  the  time,  g  is  the  constant 

gravitational  acceleration,  P  is  the  constant  density,  p  is  the  pressure, 

*  * 

and  z  =*  n  is  the  equation  for  the  free  surface.  To  derive  the  K-dV 
equation  for  the  wave  amplitude,  we  make  the  following  assumptions.  The 
channel  bottom  varies  slowly  in  the  longitudinal  direction  and  the  magnitude 
of  the  transverse  velocities  is  much  smaller  than  that  of  the  longitudinal 
velocity.  Needless  to  say,  within  the  framework  of  long  wave  approximation 
different  scalings  used  may  give  rise  to  different  equations.  However,  in  the 
following  we  shall  be  only  concerned  with  the  derivation  of  the  K-dV  equation. 
Based  upon  the  assumptions  of  the  slow  variation  of  the  channel  bottom  and 
different  orders  of  magnitude  of  the  velocities  in  different  directions,  we 
introduce  the  nondimensional  variables: 
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Figure  1 


A  cross  section  of  the  channel 
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-3/2  *  Vo  .  -3/2  *  *  * 

t  -  0  t  /(H/g)  '2  (x,y,z)  ■  (B  x  /H,  y  /H,  a  /H), 

*  *  * 

n  «  n  /h,  h  «  h  /H,  p  *  p  /(PgHK 

(u,v,w)  -  (u  /(gH)  ^2  ,  3  ^2v  /  ( gH)  ^2  ,  3  w  /(gH)  ^2  ), 
e3/2  -  l/h  »  i, 

where  L,H  are  respectively  the  longitudinal  and  transverse  length  scales. 

In  terms  of  them,  (1)  to  (7)  become 

u+3(v+w)=0,  (8) 

x  y  z 

u+uu+p+3(vu+wu)=0,  (9) 

t  x  x  y  z 

2 

v+uv+3(w+wv)  +  3p=0,  (10) 

t  x  y  z  y 

w^+uw  +3(vw  +  ww  )  +  B2  (p  +1  )  =  0 ,  (11) 

t  X  y  z  c  Z 

subject  to  the  boundary  conditions 

n  +  un  +  3(vn  -  w)  «  o  >  (12) 

t  x  y  [ 

>  at  z  =  n, 

p  =  0  J  (13) 

uh  +  3 ( vh  +  wh  )  ~  0  at  h  =  0.  (14) 

x  y  z 

We  assume  that  u,v,w,p  and  rj,  as  functions  of  t,x,y  and  z,  also 
depends  explicitly  upon  a  new  variable 
£  *  3  S(t,x) 

where  S  will  be  called  a  phase  function,  and  that  they  possess  an  asymptotic 
expansion  of  the  form 

+  (€,t,x,y ,z,3)  ~  4>0  +  +  B~24>2  +  •••  •  (15) 

Substitution  of  (15)  in  (8)  to  (14)  will  yield  a  sequence  of  equations  and 
boundary  conditions  for  the  successive  approximations  by  equating  to  zero  the 
coefficients  of  the  like  powers  of  3*  The  solution  for  the  zeroth 
approximation  is  assumed  to  be  given,  and  for  simplicity  we  assume  the 
following  one: 

<u0'v0'V  =  °'  P0  =  h0=0.  (16) 
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The  equations  for  the  first  approximation  are 


kuU  +  V1y  +  wlz  =  °' 


-  +  kp.,5  *  0, 

Ply  =  P1Z  “  0, 

subject  to  the  boundary  conditions 


+  W1  =  0 
Pi  =  ni 

v1hy  +  w1hz  =  0 


at  z  =  0 

at  h  =  0, 


(17) 

(18) 

(19) 

(20) 
(21) 
(22) 


where  k  =  Sv,  w  -  -S.  •  It  is  obtained  from  (19)  that  p.  is  a  function  of 
x  t  * 

£/t  and  x  only.  We  express 

p  «  A(£,t,x),  (23) 


and  from  (17)  and  (18),  it  follows  that 


v  +  w4  =  -Ark  /u>. 
1y  1z  5 


(24) 


Upon  integrating  both  sides  of  (24)  over  a  cross  section  D  of  the  fluid 
domain  for  constant  x,  making  use  of  divergence  theorem  and  (20)  to  (22),  we 
obtain 

(k2a(x)  -  <*i2b(x)]  «  0,  (25) 

where  a(x)  is  the  area  of  D  and  b(x)  is  the  width  of  D  in  the  plane 
z  *  0#  Assume  that  A^  is  not  identically  zero,  then  (25)  implies 

w  =  kG  ( x )  f  G(x)  =  ±[a(x)/b(x)l  ,  (26) 

where  a(x)/b(x)  is  the  mean  depth  of  the  cross  section  D.  Equation  (26)  may 
be  solved  by  means  of  the  method  of  characteristics  (Courant  and  Hilbert, 

1962)  and  the  corresponding  characteristic  equations  are 


dt/do  =*  p ,  dx/dcf  *  uG(x),  dk/d<7  =  -kUG'(x) 


dw/dO  =  dS/dO  »  0, 


(27) 


where  y  is  a  proportionality  factor.  The  solutions  of  (27)  determine  a  one- 
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parameter  family  of  bicharacteristics,  called  rays, 

t  =  X  =  X(0'V' 

where  is  constant  along  each  bicharacteristic.  We  may  choose  u  =  1  so 

that  0  -  t  and  x  =  x(t#0  )•  We  also  note  that  as  seen  from  (27)  both  gj 
and  S  are  also  constant  along  a  bicharacteristic. 

The  equations  for  the  second  approximation  are 


k  U2C  *  V2y  +  W2z  +  U1x  =  °' 

(28) 

V  +  KU1U1C  +  V1U1y  +  W1U1z 

*  kp25  *  “it  *  ft*  ■  °- 

(29) 

W' 

> 

3 

II 

>< 

<N 

a 

(30  ) 

P2z  =  “  WU  ' 

(31) 

>  the  boundary  conditions 

+  k  Ulnl5  -  w2  -  wlz  n,  + 

v  -  q 

(32) 

P2  =  n2 

/  at  z  =  0 , 

(33) 

v2hy  +  w2hz  =  "  u1hx 


(34) 


at  h  =  0  . 

Differentiating  (30)  and  (31)  in  turn  with  respect  to  y  and  z  ,  adding  and 
making  use  of  (17)  and  (18),  we  obtain 


v2p. 


■*2p< 


r2  ”  * 

It  also  follows  from  (20)  to  (22),  (30)  and  (31)  that 


■2z 
h  + 


-o)2p 


1« 

h_  =  0 


We  may  define 


where 


^2yny  +  ^2z  z 


p2  =  -♦  pi55  +  v 


at  z  =  0 
at  h  =  0. 


(35) 

(36) 

(37) 

(38) 


<t>  *  4>(t,x,y,z),  a2  =  A  (?,t,x). 


By  (35)  to  (38),  $  satisfies 

v2*  -  k2  , 


<t  =  of 
z 


at  z  =  0, 


(39) 

(40) 
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<J>  h  +  $  h  =  0  at  h  =  0.  (41) 

y  y  z  z 

Equations  (39)  to  (41)  pose  a  Neumann  problem  for  <J>  ,  and  the  condition  of 
solvability  is  seen  to  be  satisfied  as  a  consequence  of  (26).  In  terms  of 
<j>  and  A,  we  may  determine  (u«|,  ,  w1 )  by  integrating  (18),  (30)  and 

(31),  and  assuming  (u^,v^#w^)  ♦  0,  p^P^  *  0  as  £  ♦  This  implies 

(u,,  vvWl)  =  a>_1(kPl,  <42> 

We  remark  here  that  other  boundary  conditions  could  be  prescribed  and  the 
final  equation  for  A  would  be  inhomogeneous. 

Now  we  are  in  a  position  to  derive  the  K-dV  equation  for  A.  From  (28) 
and  (29), 

V2y  +  W2z  "  -  <*/w)(kp25  +  +  +  w,u1z  +  «1t  +  P1x>  '  u1x* 

(43) 

Upon  integrating  both  sides  of  (43)  over  a  cross  section  D,  making  use  of 
(26),  (32)  to  (34)  and  (38)  to  (42),  we  finally  obtain 


where 


m  A  +  m  A  +  m  A  +  m  AAC  +  m  Arrc  =  0, 
Ot  lx  2  3  £  4 


mQ  =  2b ( x ) ,  m1  =  2  a(x)/G(x), 

m  =  -(G  ^x)/  h  (h2  +  h2)  ^  ds  -  G  2 ( x)G 1  ( x) a( x) , 

2  L  X  Y  Z 

m  =  -w  1  [$  (t,x,y  ,0)  -  ( t,x,y,  0 )  ]  +  3k  G  ^xJMx) 

J  y  z  y 

-1 


(44) 


(45) 


J 


m4  *  o>  //  (V<H  dy  dz, 

D 

/  h  (h2  +  h2 )  ^  ds  is  the  line  integral  along  the  boundary  h  =  0  in 

L  x  y  z 

the  cross  section  D  (Figure  1),  (hx,  hy,  hz)  is  in  the  outward  normal 
direction  to  the  boundary  h  =  0  and  y^  ~  y^  =  b(x).  The  detailed 
derivation  of  (44)  is  deferred  to  the  Appendix.  As  seen  from  (26)  and  (27), 
along  a  ray 
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d/d o  =  3/3t  +  (dx/dt  3/3x  =  3/3t  +  G(x)  3/3x 


3,  Special  cases. 


In  this  section  we  shall  explicitly  derive  the  K-dV  equations  for 
rectangular  and  triangular  channels  with  variable  cross  section,  and  show  that 
all  previous  results  are  special  cases  of  (44). 

Case  1.  Rectangular  channel  with  variable  cross  section. 

Assume  that  the  two  vertical  walls  of  the  rectangular  channel  are  given 
by  y  =  -b^(x),  y  -  b2<x)  and  the  depth  by  z  =  -d(x),  (Figure  2).  We  also 
let  b(x)  =  [b2(x)  +  bj(x)]»  The  area  for  each  cross  section  is  a(x)  = 
b(x)d(x)  and  G(x)  =  ±d  ^2  (x).  To  be  definite,  hereafter  we  shall  only 


consider  the  plus  sign  for  G(x).  It  is  easy  to  obtain  from  (45)  that 
mQ  =  2b( x ) ,  =  2b(x)d  ^2  (x), 

m2  =  b'(x)d1yi(x)  +  d  '2  (x)d' (x)b(x)/2. 

To  determine  m^  and  m^,  we  need  to  find  satisfying 

v2*  =  k2  , 


*  oj  at  z  =  0 , 


*  0  at  z  =  -d(x), 


=  0  at  y  *  -b  (x),  y  *  b  (x). 

y  1  ^ 

2  2  2 

It  is  easily  seen  that  =  k  (z  +  d),  $  =  0  where  k  d  =  u)  .  Hence,  by 

z  y 

(45)  again  we  have 

m3  =  3k  d"  1/r2  (x)b(x),  m4  =  b(x)d3(x)k4/(3«» .  (48) 

Now  we  go  back  to  the  characteristic  equations  (27).  By  integrating 
1/ 

dx/dt  *  G( x)  «  d72(x),  the  equations  for  rays  are 

fX  .  .  /  int 


/X  G  (x)dC  *  t  -  tQ  , 


where  (tQ,  xQ )  is  the  initial  position  of  a  ray.  We  choose  Xq  =  0  and 
prescribe  S  =  -t  on  x  =  0,  then  since  S  -  constant  along  each  ray,  we 


have 

S  *  -t  =  -t  +  /x  G  1(n)dn. 
°  0 


(50) 


-10- 


A  cross 


Figure  2 

section  of  the  rectangular  channel 


For  this  choice  of  S,  S  is  time-like/  and  we  obtain  from  (50)  that 

-1  -  Vo 

0)  =  1,  W  =  G  (x)  =  d  /2  (x).  (51  ) 

It  follows  from  (44),  (47),  (48)  and  (51)  that 

2Afc  +  2d1/2(x)Ax  +  d'  1y2  (x)[b'(x)b"1(x)d(x)  +  d'(x)/2]A  ( 52 } 

+  3d  ' (x)AA^  +  d(x)  A^^/3  =  0. 

If  we  set  =  0  and  assume  that  the  channel  is  symmetric  with  respect  to 
the  plane  y  *  0,  then  apart  from  some  scaling  factor  (52)  reduces  to  the 
equation  by  shuto  (1974),  and  to  that  by  Kakutani  (1971)  and  Johnson  (1973) 
when  b*(x)  =  0.  In  terms  of  a  as  the  time  along  the  ray,  (52)  becomes 


-  Vo  .  _  .  .-1 


+  d- 


(x)[b’(x)  b  (x)d(x)  +  d'(x)/2]A 


+  3  d  (x)AA^  +  d(x)A^^^/3  =  0 


where  x  is  related  to  o  and  £  by  (50) 


=  -a  +  /x  d  1y2  (n)dn  . 

o 

Case  2*  Triangular  channel  with  variable  cross  section. 


The  two  sides  of  a  cross  section  D  of  the  triangular  channel  are  given 
by  z  =  -li^xjy  -  d(x),  z  *  U2(x)y  -  d(x)  where  ^(x)  =  d(x)/bi(x), 
i  =  1,2  (Figure  3).  The  area  of  D  is  b(x)d(x)/2  and  G(x)  =  d^2(x)  /2. 
We  find  from  (45)  that 


mQ  =  2b(x)  -  2(b1(x)  +  b2(x)],  m  -  /I  d  /2(x)b(x). 


m2  =*  /2  d  J/2  (x)  [b*  (x)d(x)  +  d'  (x)b(x)/2]/2. 


It  is  also  easy  to  verify  that 

♦  =  (k2/4)[y2  +  (z  +  d( x) )2 1 , 


satisfies  (39)  to  (41).  Hence,  by  (45)  and  (55) 


z 

A 


z  =  -d(x) 

Figure  3 

cross  section  of  the  triangular  channel 
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If  s 
2Aa  + 

where 


1 


Oj  -  ( 5/ 2/2 )d~  1/2  kb, 


m,  -  (/2  d"  1/2  (x)k3/16)  [bd3  +  d(b3  +  b3)/3]. 

4  12 

is  given  by  (50),  then  (46),  (51),  (54)  and  (56)  imply 
fi  d”  1/2  (x)  [b'(x)b_1(x)d(x)  +  d*(x)/2]A/2  +  5d“1(x)AAr 


(56) 


(57) 

+  [d(x)/4  +  (d-1(x)/12)(b^(x)  -  bt(x)b2(x)  +  b2(x))]A^  =  0, 
w  =  1 ,  k=/2d  ^2  and  x  is  related  to  £  *  3S  and  o  by  ( 50 ) 
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4.  Fission  of  solitons 


We  consider  a  symmetric  triangular  channel  with  a  cross  section  defined 


by 

b(x)  =  bQ,  d(x)  -  dQ,  x  <  x^ 

=  bi  "  V  x  >  V 

where  bg,  ,  dg ,  d^  are  constants  and  there  is  a  transition  zone  in  x^  < 
x  4  X2*  For  this  problem,  we  may  use  x  to  replace  0  as  a  variable  along 
a  ray.  By  (27), 

dx/d a  =  0 ( x )  =  d^2  (x)//2 


and  the  K-dV  equation  for  a  triangular  channel  (57)  may  be  expressed  as 

\  +  (  V2  )[b’(x)b_1(x)  +  d'(x)(2d(x))~1]A  +  (5/2/2 )d"3/2(x)AA5  (58) 

+  [/2/96)d”  1/2  (x)  [b2(x)d-1  (x)  +  12d(x)]A^  =  0, 

where  b^(x)  =  for  a  symmetric  channel.  We  set 

A  =  A*  b~  ^2  (X)  d~  ^  (x)  (59) 

and  it  follows  from  (58)  that 

A*  +  ( 5/2/2 )d‘7/4(x)b“  1/2  (x)A*A^  +  (/2/96)d"  1/2  (x)[b2(x)d_1(x)  (6Q) 
+  12d(x)]A^  =  0. 

Let  us  now  consider  a  solitary  wave  moving  from  x  -  -°®  toward  x  =  +00 
as  a  progressive  solution  of  (60)  for  b(x)  *  bQ,  d(x)  »  dQ.  We  would  like 
to  find  the  conditions  under  which  fission  of  solitons  may  take  place  after 
the  solitary  wave  moves  into  the  section  where  b(x)  *  bj,  d(x)  =  d-j .  The 
basic  ideas  used  in  our  approach  are  essentially  due  to  Tappert  and  Zabusky 
(1971),  Johnson  (1973)  and  Miura  et.  al.  (1974).  Assume  that  both  d(x)  and 
b(x)  are  constant,  and  we  transform  (60)  to  the  standard  form  of  the  K-dV 
equation 
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0 


(61) 


+  6UUC  + 


by  introducing 


U  -  p  A  / q,  C  =  qx/6 


(62) 


where 

P 


=  ( 5/2/2  )d"7/4(x)b"  1/2  (x),  q/6  =  (/2/96  )d"  (x)  [b2(x)d_1  (x)  +  12d(x)]. 

(63) 


A  progressive  wave  solution  for  U  is  given  by 

2  1/o 

U  *  (c/2)sech  [(c)  /2  (£  -  cC)/2].  (64) 

In  other  words,  if  we  prescribe  U  at  C  =  0  in  the  form  of  (64),  then  U 

will  be  a  solution  for  all  C.  Suppose  that  the  transition  of  the  channel 

takes  place  at  x  *  *  0,  and  we  prescribe 

*  *  2 

A  “  Ag  sech  at  x  =  0,  (65) 

and  from  (59), 

*  -1/o  —Vo  2 

A  *  AAbrt  2  d'2  sech  a£  at  x  =  0 . 

0  0  0 

Then 

U  *  (PgA  /q0)8ech2  <*£  at  C  =  0,  (66) 

where 

pQ  -  (5/2/2 )d‘7/4b‘ 1/2  ,  qQ/6  =  (/2/96)d^  1/2  (bjd"1  +  12dQ).  (67) 

In  comparing  (65)  with  (64),  A*  and  a  must  satisfy  the  following  conditions 

poAo/qo  “  c/2/  *^/2  “  °  ' 

so  that  (65)  is  the  initial  condition  for  a  progressive  wave  solution  of  (61) 
with  p0,  q0  given  by  (67).  It  follows  that 

P0A*/(2q0  O2)  -  1.  (68) 

After  the  solitary  wave  moves  into  the  section  of  the  channel  with  b(x) 

3  b<j ,  d(x)  -  d^,  we  assume  that  the  K-dV  equation  (61)  with  (63)  given  by 
Pl  =>  ( 5/2/0 )d“?/4b"  1/2  ,  q^/6  -  (/2/96)d‘ 1/2  [b2d +  12^],  (69) 
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could  be  solved  by  the  inverse  scattering  method  (Muira  et.  al.,  1974)  and 
consider  the  eigenvalue  problem 

d2*/d£2  +  (U(0,5)  +  A)*  =  0,  (70) 

where,  by  (62)  and  (65), 

U(0,£)  -  (p1AQ/q1 )sech2a£. 

It  follows  from  the  known  results  (Landau  and  Lifschitz,  1958;  Johnson,  1973? 
Muira  et.  al.,  1974)  that  U  will  consist  of  n  solitons  for  large  x  if 
PjAQ/(q  a2)  »  n( n+1 ) ,  n  = 

X  /a  *  -m* ,  m=1,2,*##,n, 

m 

and  the  amplitude  for  the  ra^  soliton  is  given  by 

sj  =*  -2X  m  *  1,2,***,n 
m  m 

and  by  (62)  and  ( > i ) . 

*  *  2 

A  3  q«  U  /p  =  2km  /n ( n+1 ) ,  m  =  1,2,***,n.  (72) 

m  n  1  u 

o  2 

By  eliminating  A^/cx  from  (68)  and  (71),  we  obtain 

P^qA^Pq  35  n(n+1)/2,  n  *  1,2, •••  (73) 

which  yields  a  relationship  among  bQ,  dQ,  b^  and  for  the  fission  of 

n  solitons  to  take  place  Making  use  of  (67)  and  (69)  and  letting  =  b^/d^ 

i  *  0,1,  we  may  express  (73)  as 

((Eq  +  12)/(e^  +  12)]  (efl/e  r  1/2  <d.,/d0f  11/4  =  n(n+1)/2,  n  =  1,2,«»*.  (74) 

In  the  following  we  consider  three  special  cases  upon  which  our  numerical 
results  will  be  based. 


(1) 


0  1 

From  (74),  we  obtain 


d1/d0  *  [n (n+1 )/2] 


-4/11 


n  -  1 # 2 , *  •  •  .  (75) 
It  is  clear  from  (75)  that  d^/dg  raust  be  less  than  one  for  the  fission  of 
n  >  1  solitons  to  take  place. 
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(2) 


^  =  dr 
(74)  now  yields 

[(b0/dp)2  +  \2)[h^/dQ)2  +  12]"1(b1/b0)"  1/2  =■  n(n+1  )/2  ,  n  =  1,2,3, •••  .  (76) 

For  simplicity,  assume  bQ  “  d0  =  ^ •  T^en  from  (76),  we  obtain 

13b~  1/2  (12  +  b2)”1  =  n(n+1  )/2,  n  =  1,2,3,  •••  .  (77) 

Since  the  left  hand  side  of  (77)  is  a  monotonically  decreasing  function  of 
b-j  and  equal  to  one  as  =  1,  all  the  roots  of  (77)  must  be  less  than  one 
as  n  >  1-  Hence,  the  fission  of  solitons  occurs  when  b^  <  1* 

(3)  bQ  =  br 

We  have  from  (74)  again  that 

[(d0/b0)-2  +  12] [ (d^/bQ )~2  +  12]"1(d1/d0)"9/4  =  n(n+1)/2. 

Assume  =  dQ  *  1 •  Then  it  follows  from  the  above  equation  that 

13d^1/4(12  +  d^)_1  =  n(n+1)/2,  n  =  1,2,3,  •••  .  (78) 

Equation  (78)  also  shows  that  d^  <  1  is  a  necessary  condition  for  fission  to 
take  place* 

There  are  now  several  numerical  methods  available  for  the  solution  of  the 
K-dV  equation.  Since  the  physical  models  described  by  the  K-dV  equation 
represent  situations  requiring  computations  for  large  time,  any  numerical 
method  proposed  must  meet  at  least  two  requirements.  First,  the  method  must 
yield  sufficiently  accurate  wave  amplitudes  for  many  time  steps  in  the  process 
of  computation.  Secondly,  since  the  position  of  a  wave  front  is  just  as 
important  as  the  wave  amplitude,  the  proposed  method  must  be  capable  of 
predict. in*.  Hie  position  of  a  wave  front  with  minimal  error.  The  first 
requirement  is  easy  to  meet,  if  the  method  is  conservative  (Richtmeyer  and 
Morton,  1967).  However,  it  is  much  more  difficult  to  meet  the  second 
requirement.  In  order  to  compare  results,  we  shall  use  the  partially 
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corrected  second-order  Adams-Bashforth  scheme  as  well  as  the  Hopscotch  scheme 


to  solve  (60)  numerically  for  the  three  special  cases  considered  before,  and 
we  use  n  =  3  for  all  cases* 


m  c0  '  S  • 

The  depth  of  the  triangular  channel  is  defined  by 
d0  '  x  <  X1  =  0' 

d(x)  =  ^(2/3)(d1-d0)(x2)"3(x)3  +  (l/3)(d1-d0)(x2)~2(x)2  +  dQ,  x  <  x  <  x2. 


X  >  X„, 


1  '  ”  "2* 

where  we  choose  bQ  =  =  1 '  x2  =  0.01*  From  (75 ),  we  have  =  0*5212. 


For  x  <  x^ ,  (60)  becomes 


Ax  +  ( 5/2/2 )A  A^  +  (1 3/2/96 )A^^^  =  0. 


(79) 


We  prescribe,  at  x  =  0, 


A  -  Aq  sech 


Then  by  (67)  and  (68),  we  have 


a!  =  2qfta2/p.  =  13  a2/20. 


(81) 


(82) 


‘0  "^O  '*0 

It  follows  from  (72)  that  the  amplitudes  of  the  three  solitons  are  given  by 

A*  =  13m2a2/120,  m  =  1,2,3. 
m 

*  2 

If  we  choose  A-  =  1 ,  then  <*  =  20/13,  and 

A*  =  1/6,  A*  =  2/3,  A*  =  3/2. 

The  numerical  results  are  shown  in  Figure  4. 


(2) 


d0  “  d1* 


We  define  the  width  of  the  triangular  channel  by 


b(x) 


bQ  ,  x  <  x1  *  0, 
^(2/3)(b1-b0)(x2)“3(x)3  +  ( 1/3) (b1 
l  b1  ,  x  >  x2. 


,_b0)(X2)  2(x)2  +  bo»  x1<  x  *  x2' 
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where  we  also  choose  bQ  =  dg  =  1 #  and  X2  =  0.01.  We  obtain  from  (77) 

that  b.j  *  0.0326.  The  K-dV  equation  for  x  <  0  is  still  (79).  If  we  prescribe 

*  2 

(81)  with  Ag  =  1,  a  =  20/13  again#  the  amplitudes  of  the  solitons  are 
still  given  by  (83)  as  shown  in  Figure  5. 

(3)  bQ  =  br 

The  depth  function  d(x)  is  the  same  as  in  (1)  with  =  dQ  *  1 •  From 

(78)  d^  -  0.3824.  The  amplitudes  of  the  solitons  remain  the  same#  and  the 

*  2 

results  are  shown  in  Figure  6  if  we  choose  A  =  1#  a  *  20/13. 

We  note  that  the  approximate  method  used  here  although  confirmed  by  the 
numerical  results#  is  still  formal  and  its  accuracy  may  depend  upon  the  size 
of  the  transition  zone.  In  a  subsequent  study  we  shall  justify  this  method 
and  show  that#  if  |Xj  -  x.j  I  is  sufficiently  small,  the  solution  obtained  by 
the  inverse  scattering  method  is  indeed  an  asymptotic  approximation  to  the 
exact  solution. 
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§5.  Discussion. 

In  the  derivation  of  the  K-dV  equation,  we  transformed  the  t,x- 

coordinates  to  the  a, 0^ -coordinates.  We  choose  o  =  t,  o^  *  S  since  S  is 

constant  along  a  ray.  Consider  a  simple  smooth  curve  G:  tQ  =  t0(s),  xQ  - 

x0(s),  -00  <  s  <  00  in  the  t,x-plane,  where  we  may  identify  S  with  s. 

Then  by  (49),  S,  as  a  function  of  t  and  x,  is  determined  implictly  by 

/x  G-1(n)dn  =  t  -  t  (S).  (84) 

x0(S) 

The  Jacobian  of  transformation  from  the  t ,S-coordinates  to  the  t,x- 
coordinates  is 

=  3x/3s  =  G(S)[G-1(x0)  dxQ/dS  -  dtQ/dS] . 

Hence,  J  *  0  if  dxQ/dtQ  *  G(xQ)  where  we  assume  G(S)  >  0,  and  it  follows 

that  the  rays  do  not  intersect  each  other  if  C  has  no  characteristic  direction. 

Now  we  turn  to  (53)  and  (57).  If  we  prescribe  A  as  a  function  of  x  at 

o=0,  we  have  to  change  x  to  €  according  to  (50)  where  S  *  8  An 

alternative  is  to  prescribe  S  *  x  at  t  =  0.  Then  by  (84),  S  is 

determined  implicitly  by 

/x  G_1(n)dn  =  t. 
s 

In  this  case  S  is  space-like  and 

w  =*  G ( S ) ,  k  -  G(S)/G(x).  (85) 

By  (44),  (48)  and  (85),  we  have 

2A_  +  d"  1/2  (x) [b‘ (x)b-1(x)d(x)  +  d*(x)/2]A  +  3d  ^  (S)d_1(x)AAr 
o  t 

5/2 

+  d(x)d  '  (S)A££^/3  =  0, 

for  a  rectangular  channel,  and  by  (44),  (54),  (56)  and  (8  ),  we  obtain 

2A„  +  /2d~  ^  (x)[b'(x)b~^(x)d(x)  +  d'(x)/2]A/2 

o 

+  5  d-1(x)d1/2  (S)  AA^  +  d3/2(S)  [d-1  (x)b2(x)  +  12d(x)  ]  A^/48  *  0, 
for  a  triangular  channel.  In  both  cases,  the  initial  data  can  be  explicitly 

expressed  in  terms  of  3 
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In  this  work,  we  only  consider  the  drivation  of  the  K-dV  equations.  If 


we  linearized  the  governing  equations,  and  carried  out  the  ray  method,  we 
would  get  the  same  Hamilton-Jacobian  equation  (26)  and  the  following  equation 
would  hold: 

m0At  +  m1Ax  +  m2A  =0,  (86) 

where  mQ,m^  and  m2  are  given  in  (45).  Along  each  ray,  we  may  use  x  again 
as  a  variable  and  (8  )  assumes  the  form. 

2b(x)G(x)Ax  +  m2A  =  0  (87) 

by  (27)  where  dx/dt  =  G(x).  It  follows  from  (8  )  that 
^{a  exp  /Xm2  (£)  (2b(£)G(?) )  =  0, 

along  a  ray,  and 

A  =  A0  exp[-/X  m2(C)(2b(?)G(5))'1dC],  (88) 

xo 

where  AQ  is  the  value  of  A  at  x  =  xQ,  and 

m  (£)  =  -G  (5)  /  h-r(C,y,z)(h^(C,y,z)  +  h2(C,y,z))  ^2ds 
1  L  *  y  Z 

-  g"2(Og’  (C)a(5). 

Note  that  the  attenuation  factor  in  (88)  only  depends  upon  the  geometry  of  the 
channel  and  may  be  considered  as  a  generalization  of  the  so-called  "Green’s 
law"  for  amplitude  decay  (Lamb  1932).  For  rectangular  and  triangular 
channels,  we  obtain  from  (47),  (51  ),  (54)  and  G(x)  -  d1/2(x)//2  for  a 
triangular  channel  that  the  same  Green's  law  holds: 

A  =  A0  b“  1y2  (x)d~  %  (x), 
which  is  essentially  (59). 
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APPENDIX 


We  integrate  both  sides  of  (43)  over  a  cross  section  D  (Figure  1).  By 
using  divergence  theorem,  (32)  and  (34),  we  have 


II 


( v  +  w_  ) dydz  = 

2y  2z 


h  (h2+h  )“  1/2 
x  y  z 


ds 


♦  / 
y. 


W2lz=0dy 


2^2  -  V; 


=  _U1  l  hx<Vhzr  2ds  +  /  t-"P2C  +  kUinU  "  W1zni  +  nit]z=Ody 
F  y1 

=  -//  (k/w)(kp2£  +  ki^u^  +  u1t  +  p1k)  +  ulx]dtdz 


where  yly  =  u1z  =  0.  Now  we  make  use  of  (21),  (23),  (26)  and  (38)  to  (42)  to 
compute  the  coefficients  mg  to  m4  in  (44). 
mg  =  b(x)  +  k2a(x)/w2  =  2b(x), 


m1  =  2(k /o>)  a(x)  =  2a(x)/G(x), 

m  =  -(k/to)  /  h  (h2  +  h2)“1/zds  +  (k/u) (k/u)  a(x)  +  (k/u)  a(x) 

*  _  X  V  Z  t  X 

T  * 

m  -G  *(x)  /  h  (h2  +  h2)  ds  -  G  2( x)G* (x) a(x) , 
p  x  y  z 


where  (k/w)t  ”  [G  1{x)]fc  *  0  , 

m  =  (k2/u)b(x)  +  u»  ^  /  2  $  (t,x,y,0)dy  +  (k4/o»3)a(x) 

J  zz 

y1 

=  2kG  1(x)b(x)  +  U)  1  f*2[k2  -  ♦  (t,x,y,0)]dy 

yl 

-  3kG"1(x)b(x)  -  if"1!*  <t,x,y,,0)  -  *  <t,x,y,,0)]. 

y  z  y  ' 

y 

n>4  *  <*>  /  2$(t,x,y,0)dy  -  (k2/W)  //  $dydz 

y1  D 

y 

»  u>  /  2  $(t,x,y,0)dy  -  to  1  //  ♦  ^2^dydz 

y1  D 

y  y 

-  0)  /  2  <fr(t,x,y,0)dy  -  (i)  2W2p(  t,x,y,0  )dy  +  u>  *//  (V$)2dydz 

yl  yi  D 

*  w  1  //  ( )2dydz  • 
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